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Abstract 

We analyze the frictionless motion of a point-like particle that slides under gravity 
on an inverted conical surface. This motion is studied for arbitrary initial conditions 
and a general relation, valid within 13%, between the periods of radial and angular 
oscillations, that holds for both small and high energy trajectories, is obtained. This 
relation allows us to identify the closed orbits of the system. The virtues of this 
model to illustrate pedagogically, how in a physical system, the energy is transferred 
between different modes, are also emphasized. Two easy identification criteria for 
this type of motion with potential interest in industrial design are obtained. 
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1 Introduction 



The conical shape is ubiquitous in nature. Many important large-scale natural for- 
mations such as the profile of volcanoes |jl|, 0] or the borders of tornadoes present this 
universal form ||^. On a smaller scale, different devices such as silos, loudspeakers, 
funnels, pipes, particle precipitators or centrifugal separators have been designed by 
man using this geometry to perform diverse technical operations 0, |^, ^ ^ . Don't 
let us forget other basic devices with conical design such as wash-hand basins and 
lavatories. Vortices in fluids appear in some cases as evolving cone-like structures [§], 
and recently the avalanches of two-dimensional automata modeling the surface of a 
sandpile have received considerable attention from the scientific community ||^. Al- 
though the physical mechanisms behind all these complex phenomena are different, 
the coincidence and the ubiquity of the conical form in all of them is remarkable. 

Bodies falling over inclined surfaces under the influence of gravity are also en- 
countered in every day life. Think, for instance, of a tobbogan, a roller coaster or a 
snow-board car. These systems are modeled in general physics courses as point-like 
particles sliding over those surfaces under the action of gravitational and surface 
reaction forces. 

In this work, we consider the mechanical problem of a point-like particle of mass 
m sliding on the inside of a smooth cone of semi- vertical angle (pQ, whose axis points 
vertically upward (see Fig. 1), as a first step to the understanding of those more 
complicated hydrodynamical and granular systems spiraling on the inside of a conic 
surface. This problem is proposed by T.W.B. Kibble in his well known Classical 
Mechanics textbook There, the student is asked: (1) to find the Hamiltonian 
function, using the distance r from the vertex and the azimuth angle 6 as generalized 
coordinates ; (2) to show that circular motion is possible for any value of r, and to 



determine the corresponding angular velocity luq ; (3) to deduce the angle 0o where 
the frequency of the small radial oscillations about the circular motion is also cuq. 



Setting the coordinate frame at the cone vertex, denoting the vertical direction 
as the OZ axis and using r and 9 as the generalized coordinates, we deduce that 
the Hamiltonian function for this problem is 

where Pr and po stand for the canonically conjugated momenta of r and 4>o respec- 
tively. Besides, we have assumed that the level of reference for the potential energy 
is z — 0, and the constant value of gravity is denoted by g. 

As the coordinate 9 does not explicitly appear in Tt, 9 is a. hidden variable and from 
the equations of motion we deduce that pg is a constant of the motion. Physically 
this is clear because there are only two forces acting on the particle, one being grav- 
ity and the other the force of reaction exerted by the surface which is perpendicular 
to the surface. Thus, both forces are located in the vertical plane formed by the 
OZ axis and the position vector of the particle, and therefore the momentum of 
these forces is always horizontal and hence the vertical component of the angular 
momentum is conserved. In other words pe is a constant. 

In consequence, we define the effective radial potential energy U (r) for this problem 
as 

^ ('^) ^ n 2 • 2 , + cos 00, (2) 

zmr^ sm 0o 

where L2 is the constant value assumed for pg. 

U{r) for Lz — 1 and 0o = f , is plotted in Fig. 2. Thus U{r) is a confining potential 
and the particle, no matter what its energy is, cannot escape from the cone. This 
confinement is due to the effect of gravity for large r and to the reaction of the 
conical surface as a consequence of the acquired rotational kinetic energy for small 
r. In both cases, the conservation of pg limits the transference between gravitational 



energy and kinetic energy as will be discussed in detail below. From Eq. 
deduce the radial position tq where U has its minimum Uq, 



To 



m?g siv? 0o cos 0o 



we 



(3) 



Thus, for any pair 0o and L^, there exists a unique circular motion at the distance 
r = ro, with angular velocity 



mrQ sin^ 0o 



(4) 



Now, imparting a small L^-conserving energy perturbation A to this circular motion, 
up to an energy = t/o + A, one generates an orbit with small radial oscillation 
around ro (see Fig. 2). The frequency uJj. of this oscillation is 
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Thus the radial time period Tr = — coincides with the period of the circular motion 



when Tr = —. This condition leads to 



sin-^(l/v^) ^ 35.3°, 



(6) 



which completes the answer to the three academic questions proposed at the begin- 
ning of this Section. 



In Section 0, we will study the equations of motion and the trajectories of this 
system with arbitrary energies and will obtain a general relation between the periods 
of radial and angular oscillations. In Section ^, several universal geometrical and 
dynamical relations for a general trajectory are derived. The rich energy transfer 
process acting in this system is also studied. In Section |, we state the conclu- 
sions. Finally, in the Appendix ^ we will study the necessary conditions for having 
periodicity in a general orbit. 



2 Small and Large Radial Oscillations 



2.1 Equations of the system 



As said in Section |1|, the dynamics of the particle inside the cone can be expressed 
in the generahzed coordinates {r,6) (see Fig. 1). Iff is the velocity of the particle, 
its kinetic energy T is written in these coordinates as 

T = - mv"^ = - m{f'^ + r^^^ sin^ 0o)) (7) 

where the dot means the time derivative. The potential energy V of the particle in 
the gravitational field is 

V = mgz = mgr cos 0o- (8) 
Then the Lagrangian function C for this system in these particular coordinates 
appears as 

C = T — V = - ni{r'^ + r^^^ sin^ 0o) ~ f^g^ cos 0o- (9) 
The generalized momenta are 

dC , , 

Pr = ^ = mr, (10) 

dC 

Pe = — - = siv? (pQ. (11) 
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If we substitute these new coordinates {pr,Pe) in the function 

n = rpr + Ope-C = T + V, 
the Hamiltonian Ti shown in Eq. (|I|) is obtained. 

The equations of motion of the particle are determined by the Lagrange's equa- 
tions: 
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The first of tliese equations gives us tlie evolution of the particle in the radial direc- 
tion. It yields 

f — rO"^ siv? (t>o + g cos 0o = 0. (12) 

And the second puts in evidence an invariant of motion because C is independent 
on the angular variable 6: 

mr'^ 9 sin^ (po = cte. = Lz- (13) 

This dynamical constant, L^, is the vertical component of the angular momentum. 
If we substitute the value of = ^^/(mr^ sin^ 0o) in Eq- (|T^), the radial evolution 
is uncoupled from its angular dependence which only remains present through the 
constant L^, 

r-l ^ + (?cos0o = O. (14) 

This last equation corresponds to an integrable nonlinear oscillator in the radial 
direction. After integrating this motion the angular part of the dynamics is obtained 
through Eq. ([TB|). Next we analyze the behavior of the radial oscillator: first for the 
small oscillations and then for arbitrary energy perturbations. 



2.2 The radial motion 



We can see at a glance the whole picture of the possible motions in the radial coor- 



dinate if the totality of its integral curves on the plane (r, f) are known |11 



Equation (|1^) can be easily integrated since the radial and time variables are sepa- 
rated. If we define the parameters 

A = ( ^' V. 

ymsin^o/ 
B = 5- cos 00, 

and h is the constant of the first integration of Eq. (^^, we obtain 

y + Vir) = h, (15) 



with 
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(16) 
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where f/(r) is the effective radial potential energy given by Eq. (§) and Ti is the 
total energy of the system (Eq. (|l])). 

Circular motion: Depending on the initial conditions (r, 'r)f=o, the energy con- 
stant h takes different values and a different equi-energy curve is drawn by the 
system for every h (see Fig. 3). The value h = ho for which the integral curve 
degenerates in an isolated point, occurs for the only minimum of the potential V{r), 
given by the relation 



where tq was written in Eq. (|^) and represents the only existing circular orbit of 
the system for the actual values of and 0o- By substituting = mrgCJo sin^ 0o 
in the expression (^, the angular frequency uo of the circular motion given by Eq. 
(§) is found to be 



It is curious that the angle 00 which minimizes tq for a fixed is the complementary 
angle of 4>o, written in (|), that is, 0o = 54.7°. 

General radial orbit: For h < ho motion is not possible and for h > ho the phase 
plane organizes itself as a pattern of closed curves nested around the equilibrium 
point r = To (Fig. 3). This is a consequence of the confining property of V"(r), which 
means 





g cos 00 
To sin^ 00 



lim Vir 



A 



+ 00, 



lim V{r) ~ Br 

r| oo 



+ 00. 



Hence, given the initial conditions of the system, the particle inside the cone is 
oscillating in the radial coordinate between the minimum, rmin, and the maximum, 
"^max, values of its radial trajectory (Fig. 2). Taking into account that Tq = A/B 
and that V{r) reaches the same value at these extreme points, V{rmin) = V{rmax), 
we obtain a global relation of a particular trajectory: 
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(18) 



2.3 Small oscillations 



For a value hi slightly bigger than ho, the motion of the system is represented 
approximately on the phase plane (r, r) by an ellipse around the singular point 
r = tq. In this case, r^m and Vmax are of the same order of magnitude as ro and the 
relation ([l8|) becomes simplified to 

' mm' max — ' Q- K-'-'^J 

If we put r = ro + p, with p small, and we linearize the expression ([T5|) , the equation 
of a harmonic oscillator is derived: 

p + ujIp ~ 0, (20) 

where 

?>g cos 00 

r=ro " 

as was advanced in Eq. (|^). This approximation is valid, provided the energy of the 
perturbation A = m{hi — ho) is mainly stored in the radial harmonic oscillation. 
This is verified when the condition — 4-^ ~ 1 is fulfilled. 

The relation between this radial frequency and the frequency of the circular motion 
is 

^ = 3sinVo- (21) 
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If we define the respective time periods, Tr = ^ and Tq = the latter expression 



implies that 



° sin 00, 



T 

which indicates that the angular part of the motion is faster than the radial when 
(po < (pQ. In this range, the particle will perform many revolutions around the cone 
axis for each radial oscillation. 



An orbit will represent a periodic motion, and thus, a closed trajectory in the 
four- dimensional phase space {r,6,r,6), when the ratio uJr/ujQ is a rational number 
(see Appendix In the case of small radial oscillations the perturbed circular 
orbit will be periodic if there exists a pair (p, q) of integers verifying 

2 

^ = 3sin2 0o. (22) 

This means that for each pair (p, q), with p < q, there is an angle (pQ of the cone 
opening for which the perturbed circular orbits are periodic. If this last relation 
( |22| ) is not satisfied for any pair (p, q) , the dynamics is quasi-periodic and fill up the 
torus on which it develops. As an example, if p = g = 1, then 0o = 0o — 35. 3*^ (Eq. 
(P)) and the small radial-oscillating trajectories are ellipses on the cone. 

The wave shape of the perturbed circular orbit in the plane {6, r) is not sinusoidal 
as can be checked in Fig. 4. The trajectory covers a longer 6'-angular distance when 
its radial coordinate is under the value r = ro than when it is in the upper region 
of tq. a look at the first order approximation to the orbit explains this fact. If we 
write 

r{t) ~ To + pcos{uJrt), 
then, from Eq. (|I3|) , we obtain 

e{t) (^^) (^^) sm{ujrt). 



A whole oscillation in r is made when 9 runs on the interval ^ ■ [0,27r]. If the 
perturbation in the ^-coordinate is not considered, the orbit r(9) is sinusoidal and 



satisfies r{9) < tq when 6 G 



wo 



2' 2 



Actually, a simple calculation shows that 



the small perturbation forces the dynamics to stay under r = tq on an enlarged 
6'-interval, namely ^ ■ f ~ ^ + 7^ • This shape wave deformation continues in 
the same direction by increasing the energy of the particle, in such a way that the 
motion verifies r{9) < tq for the biggest part of the ^-interval ^ ■ [0, 27r], except for 
a very narrow region where r(9) completes the oscillation and reaches rmax with a 
sharp peak in the {9, r) representation. 



2.4 Large oscillations 

For h ^ Hq, large oscillations in the radial coordinate are obtained (see Fig. 3). 
The trajectories become ovoid-like reaching the maximum radial velocity Tmax for 
r = ro: Vmax = r{ro). As Vmax > Tmin, the relation (|l|) simplifies to 

Q ,^-.2 3 
^ ' min' max — ' Q- 

Finding the relation between the radial and angular frequencies requires a more 
elaborated calculation in this case. We proceed to integrate Eqs. (0) and (|T^) 
when /i — > cx), in order to obtain the asymptotic value of uir/oje- 

We consider a half-oscillation between the extreme points, r^m and r^naxi of an 
orbit. We divide this interval into two parts: {rmini^o) and {rQ^rmax)-, which are 
covered in the time intervals: (0,ti) and (^1,^2), respectively. As the variables in 
the radial equation (p!^) are separated, we write 



dt = / , = / . 23 

2{h~V{r)) Jr^^n V2hr^ - 2Br^ - A 



This elliptic integral cannot be put in elemental functions. When /i — > 00, the 



potential V{r) behaves as 

V{r)c:^Vi{r) = — for rmin < r < ro, 
V{r) ^ V2(r) = Br for tq < r < Tmax-. 

and, Tmin and r^ax are now easily calculated, 

^l(j"min) — h =^ fmin ^ 
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V2{rmax) ^ rmax ~ ^■ 

Under these assumptions an approximated value of integral (|23|) is obtained. The 
error is basically concentrated in the region (ro — Ar,ro + Ar). However, because 
velocity is maximal in rg and tends to infinity when h ^ oo, the time invested by 
the system in that region tends to zero and thus the error vanishes. 



Substituting V{r) by Vi{r) and integrating expression (p3D in the interval {rmin, ro), 
we get 

h-^^^fi^. '-^ (24) 



Making V{r) ~ V2(r) in the region {rQ,rmax), the integral (53) is approximated by 



t2-ti^ \I^Vr^ax - ro (25) 

Thus the system spends most of the oscillation time in the upper region of the 
circular orbit r = tq, and it tends to infinity when h ^ oo. 

The angular part of the motion is derived from Equation (0), 

' d(^= 2. / (26) 
m sm 00 "'0 ^ (^) 

where the integral is performed in the angular regions (0,^i) and (^1,^2) which are 

covered when time runs in the intervals (0,ti) and (^1,^2), respectively. 

Taking into account the relation 

r'^{t) ~ 2ht^ + r^i„ for < t < ti. 



2h 



and, using the value of t = ti given by expression (0), we obtain 9i when h oo: 

Oi = , • (27) 

2 sin 00 

The behavior in the second interval (^1,6*2) is quite different. Starting from the 
relation 

'2 



t-h ~ \l— [yrmax - ro - \/rmax - r j for ti < t < t2 



r(t) can be derived. As < 1, expression (^) is approximately integrated in the 



interval (^1,^2)- We get in the limit h —>■ 00: 

02-6, - -J I, (28) 
ro\ h 

which shows that for large energies the dynamics for r > rg is projected onto a tight 
peak in the {6, r) plane. Thus most of the part of the ^-coordinate is covered when 
the particle is under the circular orbit r = tq although the system spends its time 
essentially over that circular orbit (see Fig. 5). 

Summarizing, if ^ ~ ^2 is the semi-period of the radial oscillation, the frequen- 
cies Wr and ug for large energies are 

TC 271 
ti) = ~ 

62 it/ sin 00 

which implies that, when h —>■ 00, the frequency ratio is 



— = 2 sin 00- (29) 

Ug 



Comparing this last expression with the frequency ratio ( pl| ) for small radial 
oscillations, a general behavior of this quantity can be advanced: 

^ ( V^<k<2 

— = /c sin 00 with k in the range < t It- (30) 

[ ho<h<oo 

This relation will be computationally and analytically discussed in the next section. 



3 Universal Relations of the Dynamics 



3.1 Universal equations 

As has been established in the last section, the evolution of the particle sliding inside 
the cone is governed by Eq. (14) in its radial part and by Eq. (13) in its angular 



one. These equations are expressed in a form which depends on the gravitational 
field and surface geometrical properties and, evidently, on the initial value of the 
global dynamical constants. We proceed now to write them in a universal form 
by the rescaling of radial, angular and time variables. In short, all information on 
the dynamics will be contained in these reduced equations which are independent 
of any 'external' information. To obtain the real trajectory from these equations 
it will be enough to undo the change of variables with the 'external' information: 
mass, geometry, field and initial conditions. 

If we substitute the value of the vertical angular momentum, L^, given by Eq. 
(1^) and we perform the change of variable: f = r/rg, radial and angular equations 
are reduced to 



f + 4^(l-f-3) = 0, (31) 



where uj'^^ho ^^^^ values given by Eqs. (^) and (^, respectively. 

Rescaling the time variable hj i = and the angular variable by = 6'sin0O) 

the non-dimensional equations of the motion are found: 

f + (l-f-3) = 0, (33) 
_ 0. (34) 

Remarkably, in this form, the equations apparently lose any characteristic of the 
system. Every possible trajectory on the conical surface is projected into a solution 



of these equations and, in that sense, we say that they are universaL They contain 
all the information about the dynamical behavior of the particle. 
In particular, f = 1 is the singularity representing the circular orbit and any other 
orbit runs between the extreme values, fmm and Vmax, which satisfy: < Vmin < 1 
and 1 < Tmax < oo. Expression ([T8| ) is now rewritten as 

Q ^2 ~2 

^ ' nrrtinrt ' fi 



min' max 



1. (35) 



3.2 Universal frequency relation 

From the universal equations of the dynamics, we remake now all the calculations 
presented in Subsection to show the universal relation (^). For simplicity in 
the notation, we rename the variables: r ^ r and t ^ t. 

The energy conservation becomes in the new coordinates 

^ + V(r) = E (36) 
where Vir) = r + ttt and E = — ^^—j-- Hence, in this representation, the dynamics 

' 2r^ gro cos (po ' i j 

settles in the circular orbit when E = 3/2, and, small or large oscillations are 
obtained when the value of the normalized energy E runs on the interval (|, oo). 
For an arbitrary energy G (|, oo), r^m is deduced from the equality: 

The radial coordinate as function of the time, r{t), is calculated after performing 
the integral: 

/■* , r rdr , ^, 

dt= r = r(t,E). 

For t = with Tr the semi-period of the radial oscillation, Vmax is reached. Recall 
that TmaxiE) is obtained from Expression (^). Thus, we get Tr{E) by solving the 



equation 



The angle 62 covered by the particle during a radial semi-period is given by 

rTr/2 dt 



r02 ^ 
62= d9 
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O2 = 02(E). 
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Therefore the frequency ratio obeys the following relation: 



TT 
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— sin 00 = /^(-E) sin 



where 



1 , Tr(E) 
1 2 



dt 



_TT Jo r^{t,E) 
As we calculated in the previous section, 



(37) 



(3^ 



k E 



k{E = 00) 



and as can be deduced from numerical calculations of k{E) (see Fig. 6), we claim 
that 

3 

^3 < k{E) < 2 when - < ^ < 00, (39) 
which allows us to approximate integral (|38D , valid within 13%, by 

dt TT 



/o r2(t) 2'' 
for any trajectory r{t) of arbitrary energy E. 



3.3 Universal dynamical relation 



We proceed to convert the almost universal frequency ratio (^) in a dynamical 
relation. If N represents the number of revolutions that the particle performs around 



the cone in the half of a radial oscillation, then 62 = 2ttN, and expression (pT]) is 
now written as follows: 

— = % = 4? = ^sin0o. (40) 
Taking the good approximation /c ~ 2, we obtain the universal dynamical relation 

4A^sin0o ~ 1. (41) 
In the case of small cone opening, 0o ^ 1, this simplifies to 

4iV0o ^ 1, (42) 

which establishes a simple experimental test to find out if a particle describing a 
trajectory within a conical-like surface is truly governed by the equations under 
study. That is, if the particle's dynamics is equivalent to sliding without friction on 
the surface. 



3.4 Universal geometrical relation 



Following the same line of reasoning as in Section |3.3| and motivated by the func- 
tioning of the industrial 'cyclone' gas-particle separators ||, suppose a motion 
of the system developing on a truncated cone such as that drawn in Fig. 7. The 
geometrical lengths of the orbit are the maximum and minimum radii, Rmin and 
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in the relation (jlS]), we find a universal geometrical requirement for an arbitrary 
trajectory 

^max ^ "^down taU 0o ^^^^ 

Rmin ~\~ Rmax '^9 

where Vdown is the linear velocity of the particle at the bottom of its orbit. 
Conservation of Lz implies that 

Rmax ^down 

B ~ ~v ' 



(44) 

with Vup the linear velocity at the top of the motion. Expression (^) can also been 
written as follows: 

Rmin _ < tan 00 ^^^^ 



Rmin ~l~ Rmax '^9 

These universal geometrical relations can be also used to decide if a physical phe- 
nomenon has an underlying dynamics such as that studied here. 



3.5 Energy transformation 



The kinetic energy % of the angular part of the motion can be written as 

% = - Lz 0, 

and the conservation of energy at a point P(r, 6) of a trajectory is the result of the 
balance between kinetic and potential energy, 

- ^ mr^ + ^ i^down -0) = mgAz, (46) 

where f and 9 are the radial and angular velocities at P, Az = (r — r^m) cos0o is 
the vertical distance from P to the bottom of the trajectory and Odown is the angular 
velocity at the bottom of the motion. 

If we rename = f and Vdown = Rmin^^down, Az is obtained from the latter expres- 
sion, 

?;2 7,2 / /R ■ ~ 

"down / 1 / -^^mm 



where = r sin 0o- 



p \ 2 

Umax 



The vertical distance, H — Az{Rmax), covered by the particle in its motion 
verifies 

Observe that not all the kinetic energy can be transferred to potential energy at the 
top of the orbit. A percentage of the total energy, given by the factor 7 : 
remains as angular kinetic energy in order to conserve L^. This transference is more 
efficient when the geometry of the system approaches that of a plane system because 
can be conserved by the large radial component of the trajectory This happens 
when 00 tends to |. On the contrary, the transfer of energy is nearly forbidden and 
is totally inefficient when the system tends to the cylindrical form, that is, when 0o 
goes to zero. 

The interesting energy-transfer process acting in this system, between the an- 
gular kinetic term, the radial kinetic term and the potential energy term is not 
habitual in simple mechanical systems of only two degrees of freedom. This type of 
idea typically appears later when dealing with considerably more complex systems 
of three or more degrees of freedom, such as the symmetric top. 

Note also that in other problems of a point-like particle sliding under gravity 
on surfaces such as planes or cylinders, this rich energy-transfer does not occur. 
There the transfer is as in a bullet motion, i.e., between the 'vertical' kinetic term 
and the potential energy term. In the case of sliding on the inside of a spherical 
surface the energy transfer process occurs as in the cone but with two important 
differences. First, in the sphere, no matter the value of the energy, the spatial range 
of the orbits is limited to the size of the system. And, second, when the particle 
crosses the equatorial line, the vertical reaction of the surface points downward and 



the particle may lose contact with the surface. 

For all these reasons, we believe that the cone problem has a great academic potential 
and deserves adequate attention in General Physics courses. 



4 Conclusions 



The understanding of the physical mechanisms of many complicated phenomena 
requires a drastic simplification in a first approach. Here we have isolated two basic 
ingredients, the conical geometry and the gravitational field, present in different 
natural formations such as volcanoes or tornadoes, and in a series of technical devices 
such as pipes, particle precipitators or centrifugal separators. 

We have studied the evolution of a particle sliding on the inside of a conical 
surface, with opening angle (po, under the action of gravity. This conservative system 
has two constants of motion, the total energy, Ti, and the vertical component of 
the angular momentum, L^, which makes it integrable. The initial values of both 
constants determine the trajectory described in its four-dimensional phase space, 
which in general is quasiperiodic and it fills up a 2-torus. 

Every orbit is closed when it is projected in the radial plane (r, r), running 
between the two extreme values, and Tmax- For a quasiperiodic trajectory, the 
oscillations in the angular part are not synchronized with the radial ones. It means 
that wherever the particle is released on the inside of the cone, it will come back as 
closely as desired to the same point. Only when the relation between the frequencies 
of these two (radial and angular) uncoupled modes is a rational number, p/q, we 
will have a periodic motion. This is satisfied when | = k{E) sin0O) where k{E) is a 
parameter dependent only on the product E ~ HL^'^^^ (see eqs. ( ^OD and (^)) and 



given by the expression (|38D. The range of k{E) has been numerically calculated 
in Fig. 6 and it varies in the interval ^/S < k{E) < 2, where the value k = ^/S is 
reached for small oscillations about the circular motion and k = 2 when the energy 
of the particle tends to infinity. 

Several universal relations characterizing the dynamics of the system have been 
obtained. First, the non dimensional equations have been derived by applying dif- 
ferent length and time rescaling transformations. Independently of the initial condi- 
tions, all the information about the motion is contained in these universal equations. 
Second, the relation between the radial and angular frequencies has been shown to 
depend only on the parameter k{E). It has been reinterpreted as an universal dy- 
namical relation of the number of revolutions that the particle performs around the 
cone for each radial oscillation. Third, every trajectory develops on a truncated part 
of the cone. The characteristic lengths and extreme velocities defining the orbit sat- 
isfy a universal geometrical relation. Finally, some relations of the transfer of energy 
between the orbital, radial and potential energy have been presented. 

The application of these relations to some realistic systems can be used to test if 
the corresponding underlying dynamics is understandable as caused by the combined 
action of the gravitational field and a force perpendicular to the motion, which 
could be identified as the reaction force of the conical profile associated with the 
phenomenon under study. It is our aim to explore this direction in a future work. 



A Appendix: Periodic Trajectories of the Motion 



Periodic motions of a dynamical system are represented by closed orbits. They 
play a fundamental role in understanding the topology and orbit organization of 
the phase space, but there are no general methods to find them. Here we sketch a 
procedure that simplifies the localization of periodic orbits in the present system. 
It consists of an adequate projection on a phase space of a lower dimension, where 
the property of periodicity is conserved. 

By a conserving periodicity projection we mean that if an orbit is periodic in the 
space of projection then it is also periodic in the original higher dimensional system. 
If so, there exists a bi-univocal relation between the periodic orbits of both spaces, 
and it is enough to find them in the space of projection. Thus the problem is 
enormously simplified. 



The motion of the particle inside the cone is governed by four autonomous dif- 
ferential equations of first order, given by Hamilton's equations of the system: 

r= ^ = ^, (49) 
dpr m' 



on 



dpe mr"^ sin 0o 



(50) 



du -pI 

Pr = — TT = 2 -mg COS (po, (51) 

or mr-^ sm 0o 



0. (52) 



Hence a trajectory of the system evolves in the four- dimensional phase space with 
coordinates (r, 9,pr,p0), or, equivalently, with coordinates (r, 9, r, 9) if the change of 
variables {pr,Pe) ~^ {^,9), given by Eqs. (p!0|-p!l|), is performed. 
As the dynamical equations fulfil the regular conditions of continuity and smooth- 
ness, the field is uniquely determined on each point, and the flow has no singularities 
in the region of phase space with physical significance. That is, those points whose 



coordinates verify: r>0,0<^^<27r and < 0o < f • 

A phase path will be closed, and then periodic, when it visits the same point twice 
and in consequence, an infinite number of different times: 

A trajectory is closed <^=^ 3T > ^ (r, 6, r, 6)t=o = {r, 9, r, 6)t=T- 



If we perform the projection {r,6,r,6) {f-,f), the integral curves on the new 
two-dimensional phase space are solutions of Eq. ( |T5| ) and they are closed orbits as 
can be seen, for instance, in Fig. 3. Therefore, for this particular projection, all 
the four-dimensional trajectories are projected on two-dimensional periodic orbits. 
This is a non conserving periodicity projection. 



By simple inspection of the problem it is obvious that periodic orbits of the 



equations system (^9| - ^ ) are those verifying {r{t),6{t)) = {r{t + T),6{t + T)). Take 
the origin of time, t = 0, when the trajectory passes through one of the extreme 
points, rmin or r^ax, of the radial coordinate. At these turning points, the radial 
velocity vanishes, r = 0. As the angular velocity around the cone 6 depends only 
on the dynamical constant Lz of the motion and on the radial coordinate r (see Eq. 
(p!3|)), it follows that (r, 9,r = 0, 9)t=Q = (r, 9,r = 0, 9)t=T, and the trajectory is also 
closed in the total four-dimensional space. 

Therefore, in order to find the periodic motions of the whole four-dimensional sys- 
tem we must study the conserving periodicity projection: {r,9,r,9) — > {r,9). The 
periodic orbits are in this case all those verifying that the ratio between its radial, 
Ur, and, angular, ue, frequencies is a rational number: 

Periodic motion ^ — = E. with g G Z. (53) 

uje q 



Note that this conclusion is well-established in literature |jT2[ , and is a consequence of 
the integrability of the Hamiltonian. The present four- dimensional system has two 
independent global constants, namely, the energy, TC, and the vertical component of 



the angular momentum, L^, which allow us to rewrite the Hamiltonian in action- 
angle coordinates. Hence the trajectories are 2-frequency quasi-periodic and, in 
general, they fill up a 2-torus. The angular velocities specifying the motion on the 
2-torus are ujr and a;^, and the condition for having a periodic orbit is given by Eq. 

dH). 
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Figure Captions 



1. Geometrical coordinates {r,6) defining tlie dynamics of the particle of mass 
m, which slides under the action of the gravitational field g on the inside of a cone 
with opening angle 0o- 

2. Plot of the effective radial potential energy, f/(r), for = 1, 0o = f and 
m=l. 



3. Pattern of solutions of the radial equation (JT^ for 0o = f • These are closed 
curves between Vmin and Vmax taking the maximal radial velocity for r = tq. (a) 
L, = 0.1, ro = 0.142 ; (b) L, = 1, tq = 0.660 ; (c) L, = 10, tq = 3.066. 

4. Small radial oscillations about the circular motion for Lz = 1: r ~ tq + Ar, 
Ar^ax ~ 0.05. (a) Uj. > Uq for 0o = f > 00? ''^o = 0.660. (b) Ur < ujq for 
00 = f < 00, To = 0.778. 

5. Large radial oscillations about the circular motion for = 1: r ~ tq + Ar, 
^rmax ~ 1-5. (a) ujr > for 00 = f > 00, ^^0 = 0.660. (b) Ur < ujQ for 0o = | < 0o, 
ro = 0.778. 

6. Numerical calculation oi k = — ^^^^^ ClS Sb function of rjrelation: (+) 
rjrelation = ^-'^ and (x) rjrelation = This curve is independent of 
and 00. In this case, the computation has been performed for = I and 0o = f. 



7. Characteristic lengths of a typical orbit developing on the truncated cone 
defined by Rmin, Rmax and H. 



